Apparatus and method for direct analysis of formation composition by magnetic resonance wireline logging

ABSTRACT

A wireline or logging while drilling device includes: a first coil configured to vary the intensity of a main magnetic field so as to produce phase modulation of signals emitted by spins in a sensitive volume; and a second coil configured to excite the spins in the sensitive volume in an inverse ratio of gyromagnetic constants of the spins; and receive signals from the spins displaced in frequency from the Larmor frequency of the spins by the phase modulation produced by the first coil.

CROSS REFERENCE TO RELATED APPLICATIONS

This application is a continuation-in-part of application Ser. No. 14/874,835, filed Oct. 5, 2015, entitled “APPARATUS AND METHOD FOR MEASURING VELOCITY AND COMPOSITION OF MATERIAL IN AND ADJACENT TO A BOREHOLE,” which claims priority to Provisional Application No. 62/060,321, filed Oct. 6, 2014, entitled “APPARATUS AND METHOD FOR MEASURING VELOCITY AND COMPOSITION OF MATERIAL IN AND ADJACENT TO A BOREHOLE,” the disclosures of all of which are incorporated by reference herein in their entirety. Furthermore, this application claims priority to Provisional Application No. 62/733,383, filed Sep. 19, 2018, entitled “APPARATUS AND METHOD FOR DIRECT ANALYSIS OF FORMATION COMPOSITION BY MAGNETIC RESONANCE WIRELINE LOGGING,” the disclosure of which is incorporated by reference herein in its entirety.

BACKGROUND OF THE INVENTION

This disclosure relates to the field of geological exploration, and in particular to apparatus and methods for continuous analysis of signals from magnetized geological formations by magnetic resonance wire line logging.

Continuous analysis of moving magnetized material by magnetic resonance creates periodic variation of the Larmor frequency of the moving magnetized material by modulating coils that periodically vary the main magnetic B₀ field strength. This periodically varies the impedance presented to the radio frequency oscillator power supply, which impedance function is a Lorentzian line. This methodology does not directly detect signals continuously emitted by the moving magnetized material.

Continuous creation and detection of signals emitted by moving magnetized material is disclosed by U.S. Pat. No. 6,452,390B1 to Wollin, and by the application of this technology to flow in a borehole environment disclosed in US 2016/0097664A1 to Wollin, both incorporated herein by reference for all purposes.

SUMMARY

The systems, methods and devices of this disclosure each have several innovative aspects, no single one of which is solely responsible for the desirable attributes disclosed herein.

In one embodiment described herein, a wireline or logging while drilling device includes: a first coil configured to vary the intensity of a main magnetic field so as to produce phase modulation of signals emitted by spins in a sensitive volume; and a second coil configured to excite the spins in the sensitive volume in an inverse ratio of gyromagnetic constants of the spins; and receive signals from the spins displaced in frequency from the Larmor frequency of the spins by the phase modulation produced by the first coil

In another embodiment described herein, a flow measurement device for measuring flow in or around a borehole of an earth formation includes: a magnet configured to generate a static solenoidal magnetic field with a field intensity that decreases in strength peripherally from the magnet; an electromagnet disposed around the magnet and configured to generate a time varying solenoidal magnetic field; and a radio frequency (RF) coil disposed around the magnet and configured to: generate an RF magnetic field transverse to the static solenoidal magnetic field to activate spins of nuclides having a same gyromagnetic ratio in a region of interest around the flow measurement device; and receive radio frequencies corresponding to the gyromagnetic ratio.

This disclosure provides for, among other things, a device including a permanent magnet and coils that can identify producible hydrocarbons with high spatial resolution. This disclosure also provides for, among other things, a method for directly determining productive levels during hydrocarbon wire line well logging. The method may include direct detection of a continuous signal from carbon 13 as a surrogate for hydrocarbon, enabled by continuous phase modulated dual resonance excitation of both hydrogen protons and carbon 13 atoms in a formation deep to borehole effects, and may involve electronic and circuitry modification of wire line logging tools. Logging speed can be adjusted to discriminate against non-producible surface-associated hydrocarbon, solid hydrocarbon, bitumin, and viscous high density hydrocarbon as those compounds have short T2 spin spin dephasing. Thus, signals received from carbon 13 atoms would be from producible hydrocarbon.

Analysis of the motion of an isolated spin, as described more completely below, yields an elliptic integral of the first kind which is not integrable in closed form. However, look-up tables may be used to provide the value of the integral for each value of the colatitude (flip angle) and each value of the periodic phase angle. The continuous reception of a signal with a known periodic phase angle phase modulates the emitted signal and provides the continuously detected continuously emitted signal with an extremely narrow bandwidth, markedly limiting Johnson-Nyquist noise, at known sideband frequencies displaced from the applied B₁ radio frequency field permitting continuous detection of the emitted signal in the presence of the applied constant frequency B₁ field by cross-correlation with those known frequencies.

The sine of the colatitude (“flip angle”) of the magnetization is the coefficient of the detectable amplitude of the emitted signal. In embodiments disclosed herein, this detectable amplitude is adjustable by the amplitude of the applied radio frequency B₁ field, and by the logging speed. The geometry of the applied radio frequency B₁ field defines the spatial resolution of the logging volume (vide infra).

The amplitude of the continuously received continuously emitted signal is a function of the quantity of magnetized material in the logging volume V (vide infra). The use of dual radio frequency coils permits detection of signal from low sensitivity low signal strength magnetized material by double resonance.

The logging tool can include a magnet configured to generate a static solenoidal magnetic field with a field intensity that decreases in strength peripherally from the magnet. The logging tool can include an electromagnet disposed around the magnet and configured to generate a similar but time varying solenoidal magnetic field. The logging tool can include one or more radio frequency (RF) coils disposed around the magnet. Dual RF coils can be configured to enable detection of signals from low sensitivity low signal strength magnetized material by double resonance.

In some implementations, the logging tool can include a controller configured to calculate the prevalence of selected nuclides in the logging volume V (vide infra) around the logging tool based on the received signal. In some implementations, the volume of interest can be a region in or around a borehole in an earth formation.

Embodiments of the method disclosed can include inserting the magnetic resonance logging tool into the borehole in the earth formation, the logging tool comprising a magnet, an electromagnet disposed around the magnet, and radio frequency (RF) coils disposed around the magnet. The method can include generating, using the magnet, a static solenoidal magnetic field, wherein the static solenoidal magnetic field has a field intensity that decreases in strength peripherally from the magnet. The method can include generating, using the electromagnet, a similar but weak time varying solenoidal magnetic field. The method can include generating, using the RF coils, an RF magnetic field transverse to the static solenoidal magnetic field and rotating at a Larmor radio frequency corresponding to the solenoidal field intensity in the logging volume V (vide infra) around the logging tool. The method can include receiving a signal induced in the RF coils by a magnetic field around the RF coils.

In some implementations, the method can include calculating the prevalence of selected nuclides in the logging volume V (vide infra) around the logging tool based on the received signals.

Details of one or more implementations of the subject matter described in this disclosure are set forth in the accompanying drawings and the description below. Other features, aspects, and advantages will become apparent from the description, the drawings and the claims. Note that the relative dimensions of the following figures may not be drawn to scale.

BRIEF DESCRIPTION OF THE DRAWINGS

FIGS. 1A and 1B show a cross-sectional side view and top view of an example embodiment of a magnetic resonance logging tool.

FIG. 2 shows a flow diagram of an example method for measuring a nuclide prevalence in or around a borehole in an earth formation.

Like reference numbers and designations in the various drawings indicate like elements.

DETAILED DESCRIPTION

The following description is directed to certain implementations for the purposes of describing the innovative aspects of this disclosure. However, a person having ordinary skill in the art will readily recognize that the teachings herein can be applied in a multitude of different ways. The described implementations may be implemented in any device, apparatus, or system that is capable of measuring nuclide composition of magnetized material.

The following disclosure pertains to measuring nuclide composition of material within a borehole or peripherally in a surrounding formation. The borehole may be present in an earth formation or a man-made structure. The borehole may contain a logging tool employing nuclear spin or electron spin magnetic resonance for measurement. Applications of the logging tool disclosed herein include, without limitation, hydrocarbon production, hydraulic fracturing, groundwater migration, contaminant diffusion, or detecting formation migration or tectonic plate shift.

Transient response of spin systems have usually been described by solutions to the Bloch equations which describes the magnetization as a function of both the applied magnetic fields and the relaxation effects on a phenomenological basis. Solutions to the Bloch equations require assumptions about the magnitude of these parameters and prescription of boundary conditions, creating a differential system applicable to a particular set of circumstances. To simplify this process further, the description below begins with the equation of motion of an isolated spin, modeled as gyroscopic precession creating a magnetic moment which interacts with moments created by applied magnetic fields, then introduces the total magnetization as a function of spin density, and finally limits the applicability of the equations so derived by the relationship between the times of the event sequence to the relaxation times. It is noted that non-producible hydrocarbons in general have much shorter relaxation times than producible hydrocarbons. Short T2, non-producible hydrocarbons are solid or highly viscous with signals that extinguish rapidly. Signals from long T2, less viscous producible hydrocarbons can be accentuated by adjusting logging speed. Formation permeability can be demonstrated directly by observing flow rather than empirically estimated.

A device (e.g. a wireline or logging while drilling device) can include at least two coils. A first coil (e.g. the coil 102 shown in FIG. 1A and FIG. 1B) can be configured to slightly vary a main magnetic field produced by a permanent magnet of the device, and a second coil (e.g. the coil 103 shown in FIG. 1A and FIG. 1B) can include at least one RF coil configured to produce two RF fields of intensity in the inverse proportion of the gyromagnetic ratios of the excited spins in the sensitive volume of interest, and to receive signals from these spins slightly displaced from the Larmor frequency of these spins by the phase modulation produced by the first coil.

Gyroscopic Precession

A rigid body free in space without any constraints can rotate permanently only about a principal axis of inertia”. If a rigid body rotates with speed ω₁ about a principal axis of inertia, and with ω₂=ω₃=0 about the other two principal axes, then the angular momentum vector (i.e. moment of momentum vector) {right arrow over (M)} has the same direction as the angular-speed vector {right arrow over (ω)}₁ (which is along the axis of rotation). Angular velocities of a rigid body about various axes in space, all intersecting in a point, can be compounded vectorially into a resultant angular speed about THE axis of rotation.

We infer from the last statement that the resultant angular speed ω₁ about THE axis of rotation taken as the resultant angular velocity Co′ can be decomposed into a vector sum of angular velocities Σ{right arrow over (ω_(n))}. Further, taking the rigid body as having mass symmetrically distributed about all axes through the center of mass yields a constant scalar moment of inertia I₀ about all axes, leading to the desired result decomposing the angular momentum vector {right arrow over (M)}:

$\overset{\rightarrow}{M} = {{I_{0}\overset{\rightarrow}{\omega}} = {\sum\limits_{n}{I_{0}\overset{\rightarrow}{\omega_{n}}}}}$

From Newton's equations, with {right arrow over (M_(G))} being the moment of external forces about the center of mass:

$\overset{\rightarrow}{M_{G}} = {{\frac{d}{dt}I_{0}\overset{\rightarrow}{\omega}} = {\sum\limits_{n}{\frac{d}{dt}I_{0}{\overset{\rightarrow}{\omega_{n}}.}}}}$

Thus the vector sum of the moments of a set of external forces Σ{right arrow over (M_(Gn))} equals the time rate of change of the total angular momentum I₀{right arrow over (ω)}, i.e:

${\overset{\rightarrow}{M}}_{G} = {{\sum\limits_{n}\overset{\rightarrow}{M_{Gn}}} = {{\frac{d}{dt}I_{0}\overset{\rightarrow}{\omega}} = {\sum\limits_{n}{\frac{d}{dt}I_{0}{\overset{\rightarrow}{\omega_{n}}.}}}}}$ Gyromagnetic Ratio

Taking the rigid body as having a symmetrical distribution of charge about the center of mass creates a magnetic moment {right arrow over (μ)} about any axis of rotation proportional to the angular velocity {right arrow over (ω)} about that axis of rotation where each element of charge dq at distance r from the axis of rotation creates an element d{right arrow over (μ)} of this magnetic moment {right arrow over (μ)} where, by definition of the magnetic moment,

${d\overset{\rightarrow}{\mu}} = {\left( {\pi\; r^{2}} \right)\left( {r\;\overset{\rightarrow}{\omega}} \right)\frac{dq}{2\pi\; r}}$

Each element of mass dw at distance r from this axis of rotation creates an element of the angular momentum (moment of momentum) dM of d{right arrow over (M)}=(r)(r{right arrow over (ω)})dw.

The ratio is assumed to be a constant:

$\frac{d\overset{\rightarrow}{\mu}}{d\overset{\rightarrow}{M}} = {{\frac{1}{2}\left( \frac{dq}{dw} \right)} \equiv \gamma}$ γ being a scalar constant termed the gyromagnetic ratio. Integrating, with boundary condition {right arrow over (μ)}=0 when {right arrow over (M)}=0 yields: {right arrow over (μ)}=γ{right arrow over (M)}. General Equation of Motion of an Isolated Spin

A magnetic moment {right arrow over (μ)} subjected to a field of magnetic induction {right arrow over (B)} will experience a mechanical moment (torque) {right arrow over (M_(G))} such that: {right arrow over (M _(G))}={right arrow over (μ)}×{right arrow over (B)}

In free space of permeability μ₀, the magnetic induction {right arrow over (B)} is proportional to the magnetic field intensity {right arrow over (H)}: {right arrow over (B)}=μ _(0{right arrow over (H)})

Equating the rate of change of the angular momentum to the applied mechanical moment (torque) yields

${\overset{\rightarrow}{\mu} \times \overset{\rightarrow}{B}} = {\frac{d}{dt}{\overset{\rightarrow}{M}.}}$

Multiplying by the gyromagnetic ratio γ and substituting {right arrow over (μ)}=γ{right arrow over (M)} and B=μ₀{right arrow over (H)} yields

${{\overset{\rightarrow}{\mu} \times \left( {\gamma\mu}_{0} \right)\overset{\rightarrow}{H}} = \frac{d}{dt}}{\overset{\rightarrow}{\mu}.}$

Defining γ′=γλ₀ yields the equation of motion of a magnetic moment (spin) subjected to a magnetic field intensity {right arrow over (H)}:

${\overset{\rightarrow}{\mu} \times y^{\prime}\overset{\rightarrow}{H}} = {\frac{d}{dt}{\overset{\rightarrow}{\mu}.}}$

Dividing by μ yields the instantaneous angular velocity of a magnetic moment (spin) {right arrow over (μ)} subjected to a time varying ambient magnetic field intensity {right arrow over (H)}:

${{{\overset{\rightarrow}{I}}_{\mu} \times y^{\prime}\overset{\rightarrow}{H}} = {\frac{d}{dt}{\overset{\rightarrow}{I}}_{\mu}}},$ which is a linear differential equation with constant coefficients, allowing superposition. Decomposing {right arrow over (H)} into H₀, h_(p), and h_(n)

Define an orthogonal coordinate system for {right arrow over (I)}_(μ), and {right arrow over (H)} as {right arrow over (z)}, {right arrow over (x)}, j{right arrow over (y)} where {right arrow over (x)}+j{right arrow over (y)} forms Gaussian planes everywhere orthogonal to {right arrow over (z)}. This allows decomposition of the ambient magnetic field intensity {right arrow over (H)} such that: {right arrow over (H)}=[H ₀ +h _(p) cos(Ωt)]{right arrow over (z)}+{right arrow over (h _(n))}e ^(jγ′H) ⁰ ^(t), where H₀ is a strong non-time variant ambient magnetic field intensity, h_(p) is co-aligned with H₀ and sinusoidally periodic at frequency Ω, and h_(n) is everywheres orthogonal to H₀ rotating in the local orthogonal {right arrow over (x)},j{right arrow over (y)} Gaussian plane at angular velocity γμ₀H₀=γB₀=ω₀, which is the Larmor frequency. Equation of Motion

Substituting and Rearranging in Differential Form d{right arrow over (I)} _(μ)=({right arrow over (I)} _(μ) ×γ′H ₀ {right arrow over (z)})dt+({right arrow over (I)} _(μ×[γ′) h _(p) cos(Ωt)]{right arrow over (z)})dt+({right arrow over (I)} _(μ)×γ′{right arrow over (h _(n))}e ^(jγ′H) ⁰ ^(t))dt since the vector product (cross product) distributes across a vector sum.

The first term creates a constant precession of {right arrow over (I)}_(μ) of angular velocity γ′H₀ about the {right arrow over (z)} axis.

The second term creates a periodic precession of {right arrow over (I)}_(μ) of peak angular velocity γ′h_(p) and temporal frequency Ω about the {right arrow over (z)} axis.

The third term represents a precession of {right arrow over (I)}_(μ) of angular velocity γ′h_(n) about an axis perpendicular to the {right arrow over (z)} axis, said axis rotating with angular velocity γ′H₀ in the {right arrow over (x)}+j{right arrow over (y)} Gaussian plane, a periodic Bloch-Siegert effect.

These three instantaneous angular velocities add vectorially to a resultant angular velocity, which when integrated over time, creates the locus of {right arrow over (I)}_(μ) in space.

If γ′h_(n)<<γ′h_(p)<<γ′H₀, the locus of the unit vector {right arrow over (I)}_(μ) describes a serpiginous line of colatitude θ and longitude φ on a unit diameter sphere, said sphere rotating with an angular velocity γ′H₀ in the local {right arrow over (x)}+j{right arrow over (y)} Gaussian plane. The rate of increase of θ is maximum at the poles (θ=0, π) and decreases at mid-latitude to 0.340 of maximum at the equator (θ=π/2). The governing equations are:

${{\overset{\rightarrow}{1}}_{u}X\;\gamma^{\prime}\overset{\rightarrow}{h_{n}^{\prime}}}\overset{\Delta}{=}{{\overset{\rightarrow}{\omega}}_{n} = {\frac{d\;\overset{\rightarrow}{\theta}}{dt} = {{{\overset{\rightarrow}{1}}_{u}\left\lbrack {{\cos^{2}\theta} + {\sin^{2}{\theta cos}^{2}\varphi}} \right\rbrack}^{\frac{1}{2}}X\;\gamma^{\prime}{\overset{\rightarrow}{h}}_{n}}}}$ ${{\overset{\rightarrow}{1}}_{u}X\;\gamma^{\prime}\overset{\rightarrow}{h_{p}}\cos\;\Omega\; t} = {{\overset{\rightarrow}{\omega}}_{p} = {\frac{d\;\overset{\rightarrow}{\varphi}}{dt} = {\overset{\rightarrow}{z}\gamma^{\prime}h_{p}\cos\;\Omega\; t}}}$ ${{\overset{\rightarrow}{\varphi} = {{\int{{\overset{\rightarrow}{\omega}}_{p}{dt}}} = {\frac{y^{\prime}}{\Omega}\sin\;\Omega\; t}}};{\frac{\gamma^{\prime}h_{p}}{\Omega} \cong 1.8}},{{{to}\mspace{14mu}{maximize}\mspace{14mu}{{J_{1}\left( {{vide}\mspace{14mu}{infra}} \right)}.h_{n}^{\prime}}}\overset{\Delta}{=}{h_{n}\left\lbrack {{\cos^{2}\theta} + {\left( {\sin^{2}\theta} \right)\left\lbrack {\cos^{2}\left( {\frac{\gamma^{\prime}h_{p}}{\Omega}\sin\;\Omega\; t} \right)} \right\rbrack}} \right\rbrack}^{\frac{1}{2}}}$ As  θ → (0, π); h_(n)^(′) → h_(n) ${\gamma^{\prime}h_{n}^{\prime}}\overset{\Delta}{=}{{\gamma^{\prime}{{h_{n}\left\lbrack {1 - {\sin^{2}{\theta\left\lbrack {1 - {\cos^{2}\left( {\frac{\gamma^{\prime}h_{p}}{\Omega}\sin\;\Omega\; t} \right)}} \right\rbrack}}} \right\rbrack}^{\frac{1}{2}}.\frac{d\;\theta}{dt}}} = {\gamma^{\prime}{h_{n}\left\lbrack {1 - {\sin^{2}{{\theta sin}^{2}\left( {\frac{\gamma\; h_{p}}{\Omega}\sin\;\Omega\; t} \right)}}} \right\rbrack}^{\frac{1}{2}}}}$ ${\int_{0}^{\tau}{\gamma^{\prime}h_{n}{dt}}} = {\int_{0}^{\theta}{\left\lbrack {1 - {\sin^{2}{{\theta sin}^{2}\left( {\frac{\gamma^{\prime}h_{p}}{\Omega}\sin\;\Omega\; t} \right)}}} \right\rbrack^{- \frac{1}{2}}d\;\theta}}$ ${\gamma^{\prime}h_{n}\tau} = {\int_{0}^{\theta}{\left\lbrack {1 - {\sin^{2}\theta\;\sin^{2}\varphi}} \right\rbrack^{- \frac{1}{2}}d\;\theta}}$ ${\gamma^{\prime}h_{n}\tau}\overset{\Delta}{=}{F\left( {\theta\backslash\varphi} \right)}$

which is an incomplete elliptic integral of the first kind with a periodic modular angle, and is not integrable.

Output Voltage in Receiver Coil

Based on the foregoing, continuous analysis of signals from magnetized materials may be achieved. The continuous reception of a signal with a periodic frequency phase modulates the continuously emitted signal and provides the continuously detected continuously emitted signal with an extremely narrow bandwidth, markedly limiting Johnson-Nyquist noise, at known sideband frequencies displaced from the applied B₁ radio frequency field permitting continuous detection of the emitted signal in the presence of the applied constant frequency B₁ field.

The sine of the colatitude (“flip angle”) of the magnetization is the coefficient of the detectable amplitude of the emitted signal. This detectable amplitude is adjustable by the amplitude of the applied radio frequency B₁ field, which detectable amplitude becomes a function of the known logging speed and of the unknown prevalence of the magnetized nuclide in the logging volume V (vide infra).

The amplitude of the continuously received continuously emitted signal is a function of the quantity of magnetized material in the logging volume V (vide infra). The use of dual radio frequency coils permits detection of signal from low signal strength low sensitivity magnetized material utilizing the Hartmann-Hahn ratio. Applying double resonance frequencies h_(n) ₁ and h_(n) ₂ to nuclides of gyromagnetic ratio γ₁ and γ₂ respectively in the same logging volume V (vide infra) during logging volume dwell time t with the condition of maximum energy transfer from h_(n) ₁ and h_(n) ₂ to the output signal (vide infra) γ₁h_(n) ₁ t=γ₂h_(n) ₂ t=π leads directly to the Hahn-Hartman ratio

$\frac{\gamma_{1}}{\gamma_{2}} = {\frac{h_{n_{2}}}{h_{n_{1}}}.}$ This is explained further as follows:

If the transverse magnetization {right arrow over (μ)}_(T) is defined as the projection of the magnetic moment {right arrow over (μ)} on the {right arrow over (x)}+j{right arrow over (y)} Gaussian plane, which plane is transverse to the {right arrow over (z)} axis, such that: {right arrow over (μ)}_(T)=({right arrow over (μ)} sin θ)e ^(jφ,) where θ is the colatitude of {right arrow over (μ)} with respect to the {right arrow over (z)} axis, and φ is the longitude taken from a zero meridian through the {right arrow over (z)} and j{right arrow over (y)} axes, the instantaneous angular velocity of {right arrow over (u)}_(T) in the {right arrow over (x)}+j{right arrow over (y)} Gaussian plane then is: {right arrow over (ω)}_(T)=(γ′h _(p) cos(Ωt)+γ′H ₀){right arrow over (z)} creating a phase incrementation φ of {right arrow over (μ)}_(T) at time t of

$\overset{\rightarrow}{\varphi} = {{\int_{0}^{t}{{\overset{\rightarrow}{\omega}}_{T}{dt}}} = {\left( {{\frac{\gamma^{\prime}h_{p}}{\Omega}{\sin\left( {\Omega\; t} \right)}} + {\gamma^{\prime}H_{0}t}} \right){\overset{\rightarrow}{z}.}}}$

A coil of N turns will subtend the rotating magnetization of {right arrow over (μ_(T))} such that:

$\mu_{c} = {{{\mu sin}\;{\theta\left( {\sin\;\varphi} \right)}} = {\left( {\mu\;\sin\;\theta} \right){{\sin\left\lbrack {{\frac{\gamma^{\prime}h_{p}}{\Omega}{\sin\left( {\Omega\; t} \right)}} + {\gamma^{\prime}H_{0}t}} \right\rbrack}.}}}$

By Faraday's law, the voltage induced in the coil is

$V_{c} = {{N\;\mu_{0}\frac{d\;\mu_{c}}{dt}} = {N\;{\mu_{0}\left( {{\left( {\mu\;\sin\;\theta} \right)\left( {{\gamma^{\prime}h_{p}\cos\;\left( {\Omega\; t} \right)} + {\gamma^{\prime}H_{0}}} \right)\left( {\cos\left\lbrack {{\frac{\gamma^{\prime}h_{p}}{\Omega}{\sin\left( {\Omega\; t} \right)}} + {\gamma^{\prime}H_{0}t}} \right\rbrack} \right)} + {{\sin\left\lbrack {\frac{\gamma^{\prime}h_{p}}{\Omega}{\sin\left( {{\Omega\; t} + {\gamma^{\prime}H_{0}t}} \right)}} \right\rbrack}\left( {\mu\;\cos\;\theta} \right)\left( \frac{d\;\theta}{dt} \right)}} \right)}}}$ $\mspace{20mu}{{{Since}\frac{d\;\theta}{dt}} \leq {\gamma^{\prime}h_{n}} ⪡ {\gamma^{\prime}h_{p}} ⪡ {\gamma^{\prime}H_{0}\text{:}}}$ $\mspace{20mu}{{V_{c} \cong {{N\left( {\mu_{0}\mu\;\sin\;\theta} \right)}\left( {\gamma^{\prime}H_{0}} \right)\left( {\cos\left\lbrack {{\frac{\gamma^{\prime}h_{p}}{\Omega}{\sin\left( {\Omega\; t} \right)}} + {\gamma^{\prime}H_{0}t}} \right\rbrack} \right)}},}$ where N is the number of turns in the coil, μ₀ is the permeability of free space, H₀ is the main magnetic field intensity, γ′ is the gyromagnetic ratio μ₀γ, h_(p) is the peak magnetic field intensity of the phase modulating field of temporal frequency Ω, and θ is the colatitude of the magnetic moment (spin) of magnetic field intensity μ.

The Fourier transform of V_(c) with respect to time is

${\left. {V_{c}} \right|_{\omega} = {\pi\; A{\sum\limits_{n = {- \infty}}^{+ \infty}\left\lbrack {{{J_{n}\left( \frac{\gamma^{\prime}h_{p}}{\Omega} \right)}\delta_{\lbrack{\omega - {({{\gamma^{\prime}H_{0}} + {n\;\Omega}})}}\rbrack}} + {{J_{n}\left( \frac{\gamma^{\prime}h_{p}}{\Omega} \right)}\delta_{\lbrack{\omega + {({{\gamma^{\prime}H_{0}} + {n\;\Omega}})}}\rbrack}}} \right\rbrack}}},\mspace{20mu}{{{where}\mspace{14mu} A} = {{N\left( {\mu_{0}\mu} \right)}\left( {\mu_{0}H_{0}} \right)\gamma\;\sin\;{\theta.}}}$ Phase Modulating Field

Three voltages are induced in the receiver coil; the first by h_(p) at a low frequency Ω, the second by h_(n) of radio frequency (RF) frequency γ′H₀=γμ₀H₀=γB₀=ω₀, and the third by the precession of the magnetic moment {right arrow over (μ)} (spin) consisting of a central frequency ω₀ with an infinite number of sidebands spaced about this central RF Larmor frequency ω₀ at frequency intervals Ω. These sidebands permit adjustment of h_(n) for the maximum output voltage in the receiver coil since they can be detected in the presence of the applied Larmor RF frequency ω₀ and the applied phase modulating low frequency Ω by rejecting these latter frequencies with circuit filters or lock-in amplifiers and/or by detection, heterodyning, and homodyne demodulation techniques employed in prior art radio receivers. The preferred first sideband voltage is maximized if the argument of the first sideband is adjusted so that:

${J_{1}\left( \frac{\gamma^{\prime}h_{p}}{\Omega} \right)} \cong {J_{1}(1.8)} \cong {0.582\left( {{Abramowitz},{{op}.\mspace{14mu}{cit}.\mspace{14mu} p.\mspace{14mu} 390}} \right)}$ Yielding

${\Omega = {\frac{\gamma^{\prime}h_{p}}{1.8} = {\frac{{\gamma\mu}_{0}h_{p}}{1.8} = {\frac{\gamma\; b_{p}}{1.8} = {\frac{2{\pi\left( {42.589 \times 10^{6}} \right)}}{1.8} \cdot b_{p}}}}}},$ then the side band frequency f_(p)≅23.6×10⁶·b_(p) where 2πf_(p)=Ω. The peak excursion of the magnetic moment (spin) from the plane containing H₀ and orthogonal to h_(n) is ±1.8 radians, or ±103 degrees.

Thus, Ω and h_(p) are so defined but are independent of the main magnetic field strength H₀ or Larmor frequency γ′H₀=γμ₀H=γB₀=ω₀.

WireLine Logging Application

If the magnetic moments (spins) {right arrow over (μ)} dwell in a space containing H₀ and h_(p) for a time sufficient to create significant magnetic field intensity (Slichter, op. cit. Ch.2.11, p. 51) and then dwell in a space additionally containing h_(n) for a dwell time such that nutation occurs through an angle θ=π, maximum energy is absorbed by the magnetic moments (spins) {right arrow over (μ)} from h_(n).

A medium (lattice) containing a distribution of magnetic moments {right arrow over (μ)}(spins) of dwell time t in a space in which these magnetic moments are subjected to both a strong static magnetic field intensity {right arrow over (H₀)} and co-aligned weak component h_(p), sinusoidally varying with a frequency Ω, will absorb energy from the magnetic moments by first order kinetics, opposed by random thermal motion, creating a magnetic field intensity {right arrow over (m)} such that {right arrow over ({circumflex over (m)})}{right arrow over ({circumflex over (X)})}({right arrow over (H)} ₀ +{right arrow over (h)} _(p) cos Ωt)(1−e ^(−t/T) ₁), (where

denotes both spatial vector and temporal phasor), {right arrow over ({circumflex over (X)})} being the complex susceptibility, T₁ the spin-lattice relaxation.

The logging volume V is the space occupied by the magnetic fields H₀, h_(p), and h_(n), as all previously defined. If the spins occupy the volume V for the time τ the volume logging rate is V/τ. After a time τ the spins in the logging volume are at colatitude (“flip” angle) θ and phase angle φ as governed by γ′h_(n)τ=F(θ†φ). Then

$\frac{V}{\tau} = {\gamma^{\prime}\frac{h_{n}}{F\left( {\theta\backslash\varphi} \right)}V}$ The detectable transverse magnetization is M=∫₀ ^(θ) m sin θ dθ, M=(1−cos θ). For maximum M,

${\frac{dM}{d\;\theta} = {{m\;\sin\;\theta} = 0}};{\theta = {{n\;{\pi.\mspace{14mu}{F\left( {\pi\backslash\varphi} \right)}}} = {{2{F\left( {{\pi/2}\backslash\varphi} \right)}} = {2K}}}}$ K is a complete elliptic integral of the first kind. φ is the modular angle. φ is periodic,

$\varphi = \frac{\gamma\; h_{p}}{\Omega}$ sin Ωt, a formulation not treated elsewhere.

$\mspace{20mu}{K\overset{\Delta}{=}{\int_{0}^{\pi/2}{\left\lbrack {1 - {\left( {\sin^{2}\varphi} \right)\left( {\sin^{2}\theta} \right)}} \right\rbrack^{{- 1}/2}d\;\theta}}}$ $\mspace{20mu}{{{{if}\mspace{14mu}\frac{d\;\theta}{dt}} ⪡ \Omega},{{{{then}\frac{d\;\theta}{d\;\varphi}} \cong {{o.\mspace{14mu}\sin^{2}}\varphi}} = {{1 - {\cos^{2}\varphi}}\overset{\Delta}{=}{{{m.\cos}\;\varphi}\overset{m}{=}{{\frac{1}{2\pi}{\int_{- \pi}^{+ \pi}{{\cos\left( {\frac{\gamma\; h_{p}}{\Omega}\sin\;\Omega\; t} \right)}{d\left( {\Omega\; t} \right)}}}} = {\frac{1}{\pi}{\int_{0}^{\pi}{{\cos\left( {\frac{\gamma\; h_{p}}{\Omega}\sin\;\Omega\; t} \right)}{d\left( {\Omega\; t} \right)}}}}}}}}}$ $\mspace{20mu}{{{{But}\mspace{14mu}{J_{n}(z)}} = {\frac{1}{\pi}{\int_{0}^{\pi}{{\cos\left( {{z\;\sin\;\varphi} - {n\;\varphi}} \right)}d\;\varphi}}}};({Bessel})}\mspace{20mu}$ $\mspace{20mu}{{{J_{0}\left( \frac{\gamma\; h_{p}}{\Omega} \right)} = {\frac{1}{\pi}{\int_{0}^{\pi}{{\cos\left( {\frac{\gamma\; h_{p}}{\Omega}\ \sin\;\Omega\; t} \right)}{d\left( {\Omega\; t} \right)}\mspace{20mu}{m = {{1 - {J_{0}^{2}(1.8)}} = {{1 - ({.340})^{2}} = {.884}}}}}}}};}$ $\mspace{20mu}{K = {{\int_{0}^{\pi/2}{\left\lbrack {1 - {{.884}\sin^{2}\theta}} \right\rbrack^{{- 1}/2}d\;\theta}} = {{2.5\mspace{20mu}{V/\tau}} = {{\left\lbrack \frac{\gamma^{\prime}}{2(2.5)} \right\rbrack{Vh}_{n}} = {{\left\lbrack \frac{\gamma}{2(2.5)} \right\rbrack{Vb}_{n}} = {{\left\lbrack \frac{2{\pi\left( {42.6 \times 10^{6}} \right)}}{2(2.5)} \right\rbrack{Vb}_{n}\mspace{14mu}{Tesla}\mspace{20mu}{V/\tau}} = {{53.5 \times 10^{6}{Vb}_{n}\mspace{14mu}{Tesla}} = {535 \times 10^{9}{Vb}_{n}\mspace{14mu}{Gauss}}}}}}}}}$ Exemplary Magnetic Resonance Well Logging Tool

With this background, a logging tool implementing these principles can be designed.

FIGS. 1A and 1B show a cross-sectional side view and top view, respectively, of an example embodiment of a magnetic resonance well logging tool 100. The tool 100 includes a magnet 101 having north and south poles, an electromagnet 102, and radio frequency (“RF”) coils 103 and 104. The tool 100 is configured to insert into a borehole 105 of an earth formation 106 suspended by a wireline and withdrawn creating relative motion between the tool and the formation.

The tool 100 can include a magnet 101. The magnet 101 can generate a static magnetic field around the flow measurement device 100 that is generally solenoidal in shape. The strength of the static magnetic field generated by magnet 101 can decrease in strength peripherally from the magnet. The magnet 101 can be a permanent magnet. The magnet 101 can be of any size or shape appropriate to the dimensions of the borehole and the measurement application for which it will be used, including a long cylindrical shape.

The logging tool 100 can include an electromagnet 102. The electromagnet 102 can be configured to generate a time varying magnetic field around the tool 100 that is generally solenoidal in shape. The electromagnet 102 can be of any type suitable for generating a solenoidal magnetic field of roughly the same relative distribution as the static magnetic field. The electromagnet 102 can take the shape of a helical coil surrounding the magnet 101. The strength of the time varying magnetic field generated by the electromagnet 102 can vary linearly with the application of a time varying current through the electromagnet 102. The strength of the time varying magnetic field generated by the electromagnet 102 can be relatively weak compared to the strength of the static magnetic field; for example, a fraction of a Gauss to several Gauss compared to several thousand Gauss in the static magnetic field. The frequency of the time varying magnetic field generated by the electromagnet 102 can be very low compared to the frequency of the RF magnetic field generated by the RF coils 103 and 104.

The tool 100 can include radio frequency (“RF”) coils 103 and 104. The RF coils can be configured to continuously generate a time varying magnetic field h_(n) transverse to the H₀ static solenoidal magnetic field generated by the magnet 101. The RF coils 103 and 104 can be a single coil or made up of multiple coils. The RF coils 103 and 104 can be a single coil configured to continuously both create and receive signals from time varying magnetic fields. Alternatively, separate RF coils 103 and 104 can be employed continuously as a transmit coil and a receiving coil, respectively. The RF coils 103 and 104 can be of a birdcage or quadrature design, or any other coil, cage, or antenna structure suitable for directing an RF magnetic field outward from the tool 100. In one embodiment, the RF coils 103 and 104 comprise a transmitting birdcage coil disposed around the magnet and configured to continuously generate the RF magnetic field, and a receiving birdcage coil disposed around the magnet and configured to continuously output the received signal. In another embodiment, the RF coils 103 and 104 could comprise a first birdcage coil disposed around the magnet and configured to operate at a first frequency, and a second birdcage coil disposed around the magnet and configured to operate at a second frequency. The RF coils 103 and 104 can be tuned to resonate at the desired frequency of operation. The RF coils 103 and 104 can have a high quality factor (“Q”) for improved efficiency of transmission and receipt of time varying magnetic fields. The RF coils 103 and 104 generate a time varying h_(n) magnetic field rotating at the Larmor frequencies of the selected nuclides (e.g. H1 and C13), ω₀=γh₀ corresponding to the H₀ static solenoidal magnetic field intensity h₀ in the logging volume V around the logging tool 100 in the ratio of the gyromagnetic constants of the selected nuclides. The logging volume V can take the shape of a segment of a barrel shaped prolate or oblate ellipsoidal surface having a finite thickness generated about the centerline of the logging tool 100. The RF coils 103 and 104 can be configured to output received signals induced by the magnetic fields about the RF coils 103 and 104.

Multiple RF coils 103 and 104 can be employed to measure the relative abundance of two different nuclides in the logging volume V; for example, detection and measurement of both ¹H and ¹³C by the same tool 100 in the same logging volume V. Exciting RF coils 103 and 104 with frequencies selected for two different nuclides permits implementation of double resonance yielding the prevalence of, e.g., ¹H and ¹³C, thereby implementing a multiphase “cut meter” or wireline logging modality. Such a system may be valuable, for example, for measuring the relative abundance of hydrocarbons versus brine in the logging volume V surrounding the tool. Such information can be valuable when searching for, e.g. hydrocarbon deposits.

In operation, the magnet 101 and electromagnet 102 can create a strong solenoidal field with a very weak slowly time varying component in the borehole 105 and in the surrounding invaded and non-invaded earth formation 106 beneath any mudcake. The RF coils 103 and 104 can create an adjustable RF magnetic field essentially orthogonal to the solenoidal field. The RF magnetic fields rotate at the Larmor frequency and causes the spins to nutate with increasing colatitude angle with respect to the strong magnetic field, permitting reception of signal induced by a magnetic field created by the nutating spins at sideband frequencies displaced from the applied Larmor frequencies. Selecting the Larmor frequency selects the depth of investigation, where the depth of investigation represents the distance from the tool 100 to the region of interest of spin activation about the logging tool. That is, for a given strength of the static solenoidal magnetic field H₀, adjusting the frequency of the RF magnetic field generated by the RF coils 103 and 104 controls at what distance from the tool 100 spins of the same gyromagnetic ratio will be affected. An RF magnetic field having a higher frequency will activate spins of the same gyromagnetic ratio relatively closer to the logging tool 100, and an RF magnetic field having a lower frequency will activate spins of the same gyromagnetic ratio relatively further from the tool 100. The strength of the slowly varying component of the solenoidal field is adjusted for optimum signal reception of the preferably first sideband frequencies from the depth of investigation.

The conditions generated by the tool 100 create both nutation, with increasing co-latitude (flip angle), and periodic phase modulation of the rotational Larmor frequency of the spins. The received signal then induced in the RF coils 103 and 104 is sinusoidal with slowly varying frequency whose Fourier transform yields a central Larmor frequency and discrete side bands displaced by the phase modulation frequency. These side bands can then be detected in the presence of the strong continuous applied Larmor excitation frequency field by means known in the art including heterodyne frequency shift, homodyne detection, and cross correlation in a lock-in amplifier. Adjusting the strength of the applied Larmor excitation frequency field B₁ for maximum received signals yield a function of the prevalence of the spins and of the logging speed of the selected nuclides at the depth of investigation in the logging volume V.

The radio frequencies supplied to each coil 103 and 104 are in the frequency ratio of the gyromagnetic ratios of the detected nuclides to excite the same logging volume. These frequencies can be chosen so that the logging volume so selected is peripheral to borehole effects. The amplitude of each supplied radio frequency is inversely proportional to the gyromagnetic ratios of the detected nuclides to permit cross excitation of the less prevalent by the more prevalent nuclide by an appropriate resonance technique (the Hahn-Hartmann effect is a non-limiting example of such a technique). The amplitudes of these supplied radio frequencies are adjusted to produce maximum signal to account for variations in logging speed. The current in electromagnet 102 can be adjusted to maximize the value of the preferably first side band about the Larmor frequencies from the received signals induced in the RF coils 103 and 104.

FIG. 2 shows a flow diagram of an example method for measuring prevalence of selected nuclides in or around a borehole in an earth formation, using the techniques described above. In block 310, a tool or device (e.g. a wireline or logging while drilling device described herein) is inserted into a borehole. In block 320, the device generates a static solenoidal magnetic field. In block 330, the device generates a time varying solenoidal magnetic field. In block 340, the device generates an RF magnetic field transverse to the solenoidal magnetic field. This field may, for example, excite spins in a sensitive volume in an inverse ratio of gyromagnetic constants (e.g. gyromagnetic ratios) of the spins. In block 350, the device may measure a received signal induced in an RF coil of the device. The received signal may be, for example, from the spins, and may be displaced in frequency from the Larmor frequency of the spins by the phase modulation produced by the electromagnet. The received signal may, for example, have a radio frequency corresponding to the gyromagnetic ratio.

The foregoing disclosure is equally applicable to nuclear and electron magnetic resonance. Furthermore, the measurement of prevalence of nuclides by the device described in this disclosure is applicable not only to liquid or gas, but to other materials, such as mixtures, slurries, aggregates, blowing particles, viscous plastics as well as to solid materials.

Preferred embodiments of the invention have now been described. It will be appreciated by those skilled in the art that such embodiments are intended to exemplify the invention. Various other embodiments of the invention will be apparent, which fall within the spirit and scope of the invention. 

What is claimed is:
 1. A wireline or logging while drilling device, comprising: a first coil configured to vary the intensity of a main magnetic field so as to produce phase modulation of signals emitted by spins in a sensitive volume; and a second coil configured to: excite the spins in the sensitive volume in an inverse ratio of gyromagnetic constants of the spins; and receive signals from the spins displaced in frequency from the Larmor frequency of the spins by the phase modulation produced by the first coil.
 2. A flow measurement device for measuring flow in or around a borehole of an earth formation, comprising: a magnet configured to generate a static solenoidal magnetic field with a field intensity that decreases in strength peripherally from the magnet; an electromagnet disposed around the magnet and configured to generate a time varying solenoidal magnetic field; and a radio frequency (RF) coil disposed around the magnet and configured to: generate an RF magnetic field transverse to the static solenoidal magnetic field to activate spins of nuclides having a same gyromagnetic ratio in a region of interest around the flow measurement device; and receive radio frequencies corresponding to the gyromagnetic ratio.
 3. The flow measurement device of claim 1, wherein the second coil is configured to excite frequencies in the frequency ratio of the gyromagnetic ratio.
 4. The flow measurement device of claim 1, wherein the second coil is configured to supply frequencies having an amplitude that is inversely proportional to the gyromagnetic ratio.
 5. The flow measurement device of claim 1, wherein the second coil is configured to excite spins of a first nuclide and a second nuclide. 